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Abstract 

In the real-time thermal field theory, the pion self-energy at finite temperature and density 
is evaluated where the different mesonic and baryonic loops are considered. The interactions of 
pion with the other mesons and baryons in the medium are governed by the effective hadronic 
Lagrangian densities whose effective strength of coupling constants have been determined from 
the experimental decay widths of the mesons and baryons. The detail branch cut structures 
of these different mesonic and baryonic loops are analyzed. The Landau cut contributions of 
different baryon and meson loops become only relevant around the pion pole and it is completely 
appeared in presence of medium. The in-medium spectral function of pion has been plotted for 
different values of temperature, baryon chemical potential as well as three momentum of the 
pion. A noticeable low mass probability in pion spectral function promise to contribute in the 
low mass dilepton enhancement via indirect modification of p self-energy for tttt loop. 

1 Introduction 

The recent reviews mm claim that the calculations based on hadronic many-body theory are very 
successful to describe the low mass dimuon enhancement measured by the NA60 collaboration [3], 
where a strong broadening of rho meson are essentially concluded as a main reason. There are 
two sources which are generally believed to be responsible for this broadening of the p meson. 
One is coming from the direct interactions of p with the other thermalized mesonic mm and 
baryonic [SEllH] constituents of the medium. Another is originating from the indirect interaction 
of p with medium via pion cloud, which may be modified through interactions with the surrounding 
medium lattniiii]. The latter source indicates that the in-medium modification of pion may create 
an effect on p meson propagation in the medium, for which the rate of low mass dileptons are 
expected to be enhanced. In this context, the present manuscript is intended to investigate the 
in-medium modification of pion propagation due to strong interaction with the other thermalized 
mesons and baryons in the medium. For this purpose the in-medium pion self-energy for different 
mesonic and baryonic loops are calculated in the real-time thermal field theory (RTF) and their 
detailed branch cut structures are explicitly analyzed. With respect to earlier investigations m 
[laiiiiiisiiisiiiTiiis], one of the major contribution of this present article is that an extensive set 
of baryonic loops are taken to estimate the in-medium self-energy of pion. The na and irp loops 
have been taken as mesonic loops because the resonances a and p have to be considered in effective 
hadronic model to reproduced the experimental data of tt-k scattering cross section m- 

The manuscript is organized as follows. The expression of thermal propagator for vr meson is 
explicitly derived in the next section and then the detail calculations of pion self-energy at finite 
temperature and density are addressed in Sec. (3). The baryon and meson loop calculations are 


1 


separately discussed in two subsection of the Sec. (3). The detailed numerical results are discussed 
in Sec. (4) and at last section the intention of the article is summarized. 


2 Pion propagator in the medium 

In the real-time formulation of thermal field theory, any two point function like propagator acquires 
2x2 matrix structure. The free propagator matrix of pion at finite temperature can be defined as 

Dab{k) = i j d"^xe*^^(rc7r(x)7r(0))a6 (1) 


where the subscripts a, 5(= 1, 2) are thermal indices, 7r(x) is pion field and Tc denotes time ordering 
with respect to a symmetrical contour in the plane of the complex time variable m- The notation 
(O) is used to represent ensemble average of any operator O, (O) = Tr^^^§-, where Tr indicates 
trace over a complete set of states. With the help of the Dyson equation, the complete in-medium 
vr-propagator can be expressed as 


Dabik) = - D^^Xk)ncd{k)Ddbik) 


(2) 


where (k) and nab(A:) are free propagator and self-energy matrices of pion at finite temperature. 
By the diagonalization procedure [2T] the thermal indices (a, b, c, d) can be removed and it gives 
the following equation : 

D{k) = (fc) - {k)Ti{k)D{k) (3) 

where the quantities with bars represent the corresponding diagonal components. The diagonal 
element of free propagator is expressed as 


D^^\k) 


-1 

k2 _ 


(4) 


which is exactly similar with the vacuum free propagator. Using Eq. ([1]) in ([3]) we get the diagonal 
element of complete pion propagator. 


D{k) 


-1 

k2 _ ^2 _ n(/c) 


(5) 


whose imaginary part provides the Breit-Wigner type structure of pion spectral function. 


A.,r{k) = lmD{k) 


—lmll{k) 

{A:2 -ml- Ren(A:)}2 + {Imn(fc)}2 ' 


( 6 ) 


Here the information of medium is embedded in the imaginary and real part of the pion self-energy 
(ImH and ReH) which will be derived in next subsection. 


3 Pion self-energy in the medium 

In this section, the different meson and baryon loop diagrams (shown in Fig. [1]) of vr meson, have 
been evaluated in order to investigate the in-medium modification of pion in hadronic matter. The 
calculations of meson and baryon loops are separately addressed below. 
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Figure 1: Self-energy diagrams of pion for meson loops (vr^^) is shown in diagram (a). The baryon 
loops (NB) of pion self-energy can be represented in two possible ways - (b) and (c). 


3.1 Meson loops 

From the lowest order xPT) th® estimated vrTr cross section in free space is well in agreement 
with the experimental data up to the center-of-mass energy ^/s = 0.5 GeV. Beyond this value 
of ^/s, the a and p resonances play an important role to explain the experimental data. By the 
unitarization technique, the a and p resonances can be generated dynamically [22] in the amplitude. 
An alternative way nais], which is followed in the present paper, is to incorporate these resonances 
by using the effective Lagrangian densities for Trvrcj and vrvrp interactions: 

C = QpP^- TT X ■ TTa, (7) 

where the coupling constants (g^ = 5.82 and Pp = 6) are fixed from their experimental decay widths 
of a and p mesons in their vrvr channels [23| . Serious defect of phenomenological interactions is that 
they often do not respect constraints from chiral symmetry. Owing to the appearance of a and 
p resonances in tttt scattering cross section, we have considered na and np loops to calculate in¬ 
medium pion self-energy where the interaction part is governed by the effective Lagrangian densities, 
mentioned in Eq. dZD- This pion self-energy function for ttcj and np loops has diagrammatically 
been represented by the FiglTJa). The 11-component of pion self energy for the 7r<l> loop (where 
denotes cr or p in generalized way) is given by 

r pi] 

nM(^) = * J - l,m^) (8) 


where the vertex factors. 


for the TTfT loop, and 


L{k,l) 



L{k,l) 



[k 


1 ) mi + k^ 


(9) 


( 10 ) 


for the Tip loop have been obtained from the effective hadronic Lagrangian, given in Eq. ([7|). The 
L)^^(/,mjr) and D^^{k — l,m^) in (|8|) are the 11-component of the scalar propagator for vr and 
meson respectively. Eor example. 




P 


—-h 2mni5{P 

771 “ -|- ip 



( 11 ) 
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where ni = —- is Bose-Einstein distribution function of pion with its on-shell energy, ui = 


P Similar to the propagator matrix, the self-energy matrix can also be diagonalized into 

a single component and the diagonal element and 11 component are related as [SUED 


ImnM(^) = tanh 




ImllMik) 


RenM(^) = . 


( 12 ) 


After doing the Iq integration of Eq. ([8]) and then using the Eq. the diagonal element of 
the pion self-energy matrix is obtained as 


^M{k) = 


dH 


1 


(1 -h ni)Li + UuLs 


+ 


—niLi -|- UuL^ 


(27r)3 Aujiuju [ko -uji-u!u + iv^^iko) ko - uji + lVu + iv^^iko) 


+ 


niL2 — UuLs 


+ 


-niL2 - (1 + nu)Li 


(13) 

ko + toi - uju + iv^^iko) ' ko + (jJi + uju + im{ko)\ 

where Lj,z = 1,..4 denote the values of L{lo) for lo = uJi,—uJi,ko — uJu,ko + uJu respectively and 
~ is Bose-Einstein distribution function for resonance with its on-shell energy, uju = 


{k — ly + m\. The imaginary part in the relevant Landau and unitary cut regions are respectively 
given below 


ImllM = — 


and 


ImllM = — 


e(go) 

167r|fc| 

e(^o) 


167r|A:| 


duJiL2{n{uii) - n{uiu = ko+ui)} 


duiLi{l + n{u!i) + n{ko - uji)} 


Si 


(14) 

(15) 


where the integration limits ujf = ^(/cqT I^I^tt), = ^(—1^1144-) with 114 ~ \j'^ ~ 
and si = k'^ — m\ -|- These two regions of branch cuts are k < ko < \lk‘^ + (m^ — 


kP' + (m$ -|- TTijr)^ < ko < oo respectively. 

The real part of the self-energy can be easily read off from (jl3p in terms of principal value 
integrals and we do not write them here. 

In realistic scenario, the a meson should not be treated as stable particle as considered in our 
previous framework as it exhibits a broad spectral function in vacuum. Hence, to take into account 
its broad width, the Eq. (fT^ has been convoluted as (see e.g. [23l [28l [8l [29]1 


IlMik,m^) = 


+ \2 


ri'm.t) 

'(m-r 


dp^ A„{p) IlM{k,p) 


/ 


+ \2 


ri'mt) 

'(my)2 


dp'^ A„{p) 


(16) 


where IiM{k,p) is the narrow-width expression self-energy given in Eq. (1131) . with m^j replaced by 
p; Afj{p^) is the spectral density: 


= -Im 


TT 


-1 


p2 _ jy^2 _|_ [p) 


ro-(p) is the vacuum spectral widths of the mesons: 

r,(p) = 


2 \ 1/2 


32np 


p^ 


(17) 


(18) 


In the integration limits, = mg- ± 2rD, with Tj) = ro-(p = 


rrin 
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3.2 Baryon loops 


Next we will calculate the baryon loops of pion self-energy in the same framework of RTF. Besides 
the meson loops in the medium, the pion propagator may also be undergone via different inter¬ 
mediate NB loops, where B stands for different higher mass baryons including nucleon itself. In 
this work we have taken an extensive set of 4-star baryon resonances with spin one-half and three- 
half. These are A^(940), A(1232), A^*(1440), A^*(1520), iV*(1535), A*(1600), A*(1620), iV*(1650), 
A*(1700), A^*(1700), A*(1710), ^"*(1720) where their mass (in MeV) are presented inside the 
brackets. The direct and cross diagrams of pion self-energy for NB loops have been represented in 
the diagram [IKb) and (c). 

The 11-component of the in-medium pion self-energy for the NB loop is given by 

n^i(/c)=i f- ak,mB) (19) 

where E^^{l,m]y) and E^^{1 — ak,mB) are respectively scalar part of the nucleon and baryon 
propagators at finite temperature. In RTF its expression is as follows : 

E^^{l,mN) = 72- ^~r~- - 2mN{lo)6{l‘^ -m%)] N{lo) = nle{lo) + n^e{-lo) 

E — + ir] 


J_ f ^ n+ _ l-rii 

2u)i \lQ-uJi + iri lo-uji- ir] Iq + uji - irj 


n, 


I 0 +UJ 1 + ir] 


( 20 ) 


where ip/j ) _|_ 7 Fermi-Dirac distribution for energy oji = \Jl‘^ + and the ± 

sign in the superscript of n; refer to nucleon and anti-nucleon respectively. Here /tat is the chemical 
potential of nucleon which is supposed to be equal with the chemical potentials of all the baryons 
considered here. The two values of a in (|19l) correspond to the direct and crossed diagrams shown 
in Fig. [1] (b) and (c) respectively which can be obtainable from one another by changing the sign 
of the external momentum k. 

Let us first discuss diagram (b) for which a = -1-1. After integrating over in Eq. (|19l) and then 
using the similar kind of relation like Eq. m, the diagonal element of the in-medium self energy 
can be expressed as 


IlBik) 


dH 


(27r)3 


(1 - nj)Li - ^ nlLi - n+Lj 

ko - LVi - Uju + iv^^iko) ko - UJI + Uju + irj^iko) 


-Hi L 2 + L 3 ^ Hi L 2 + (-1 -F n+)L4 

ko + uji - oju + irie{ko) ko+iOi + Uu + imiko) 


( 21 ) 


where are also Fermi-Dirac distribution functions for baryon and anti-baryon with uJu = 

\J {I — k)^ -|- m'^B = l,--4 denote the values of L{lo) for lo = u!i,—uJi,ko — uiu,ko + uJu 

respectively. For the diagram (c) in Fig. [H a similar kind of expression like Eq. (EH) can also be 
received just by putting a = — 1. The imaginary part of the third term of ( 12 ip will be non-zero in 

the Landau region k < ko < \Jk‘^ + {rriB — m]\r)^ and only this contribution is significant for pion 
spectral function because pion pole is situated in this region while unitary cuts are far from it. 
Adding the relevant Landau cut contributions of diagram (b) and (c), the total contribution from 
the baryon loops is given by 


Im IlBiko,k) 


_^/cn") f^i 

-^ du}i[Li{a =-l){-n+{uji) + n+{uiu = ko + uii)} 

167r|A:| 

+ L2{a = -hl){-n_(a;i) -b n_(a;„ = ko + uJi)}] 


( 22 ) 
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where iof" = ± |^|hFiv) with Wn = S% = k'^ — m\ + m^. For the B = N 

(i.e. NN loop), the diagram (b) and (c) are identical, hence only one of them has to be included. 

Similar to mesonic case, the real part of Imlls can also be extracted from Eq. (1221) by taking 
its principal value integral. The thermal part of Refls or RellM is the main interest to estimate 
the mass shift of pion, which can be reflected from the shifting of peak position in the pion spectral 
function. 


Baryons 


Ib 

Ttot 

I'b^Nti (B.R.) 

flm^ 

A(1232) 

3 + 

2 

3/2 

0.117 

0.117 (100%) 

15.7 

iV*(1440) 

1 + 

2 

1/2 

0.300 

0.195 (65%) 

2.5 

iV*(1520) 

3- 

2 

1/2 

0.115 

0.069 (60%) 

11.6 

iV*(1535) 

1 - 
2 

1/2 

0.150 

0.068 (45%) 

1.14 

A*(1600) 

3 + 

2 

3/2 

0.320 

0.054 (17%) 

3.4 

A*(1620) 

1 - 

2 

3/2 

0.140 

0.035 (25%) 

1.22 

iV*(1650) 

1 - 

2 

1/2 

0.150 

0.105 (70%) 

1.14 

A*(1700) 

3- 

2 

3/2 

0.300 

0.045 (15%) 

9.5 

A* (1700) 

3- 

2 

1/2 

0.100 

0.012 (12%) 

2.8 

A*(1710) 

1 + 

2 

1/2 

0.100 

0.012 (12%) 

0.35 

A* (1720) 

3 + 

2 

1/2 

0.250 

0.028 (11%) 

1.18 


Table 1: From the left to right colnmns, the table contain the baryons, their spin-parity quantum 
numbers J^, isospin Ib, total decay width Ttotj decay width in Ntt channels Tb->-Nw or T^ in 
Eq. (j26p (brackets displaying its Branching Ratio) and at the last coupling constants f 


To calculate L{k,l) for baryon loops we have used the effective BNn interaction Lagrangian 
densities [25] 


C 


I -^5 I + h.c. for Jb = \ , 


(23) 
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Figure 2: Imaginary (upper panel) and real (lower panel) part of pion self-energy for vrcr (solid line) 
and Tip (dashed line) loops. In the imaginary part of the pion self-energy, two separate regions of 
Landau and unitary cuts are distinctly revealed. The dotted line represents the vrcr loop contribution 
after foding by the broad vacuum spectral function of a. 


where Jb and P are spin and parity quantum numbers of corresponding baryon, B. The effective 
strength of the coupling constants f for different BNtt interactions can be determined from the 
experimental vacuum widths of corresponding B Ntt decays. From above Lagrangian densities, 
one can easily deduce 


L{k,l) = “ — PmB)^(jl + rriN)] 

\ ^TT / 

= [2(k ■ l)'^ — a(k ■ l)k‘^ — + mNiTiB)] 

Vm^/ 


for Jb = and 


L{k,l) = - 


f 


Tr 


2>m\ 


(24) 


(/ -F m7v)(/ -aJf: + PmB)k^ky 

- akYil - akY + - akY -{I- akYlY 

oTTl^ oTTIq 

[misimB + 1“^ — a{k ■ l)][{l ■ k — ak'^Y — k‘^m%] (25) 

\ ^TT / 


for Jq = 1^. Following the Ref. [T], a soft monopole form factor F = (A^ — m^)/(A^ -|- fc^) with 
A = 0.3 GeV has been multiplied in the ttNB vertex to take into account its finite size. 

The Lagrangian densities in (l23]i are not displaying its isospin structures. For ^ and 

= 1^, these isospin structures should be v/^r • Trip and ipT ■ Trip respectively, where T and 
r stand for the usual spin 3/2 transition and Pauli operator. These isospin structures provide 
appropriate isospin factors, which have to be multiplied with the expressions of corresponding NB 
loop diagrams. The isospin factor for NN or NN* loops is Y^nn,n* = 2 and for the NA or NA*, 
it is Y^naa* = ^3- 

Next we have calculated vacuum width of different baryons in the Ntt decay channel to fix their 
corresponding coupling constants f With the help of the Lagrangian densities ([7]), vacuum 
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Figure 3: Imaginary (upper panel) and real (lower panel) part of total pion self-energy coming from 
the meson loops are presented for different set of temperature T and three momentum k. 


decay width of baryons B for Nn channel can be obtained as 


Tb = 


“IJB ^ 

-^A,A*^7rjV 

2,7B + 1 


— ) + rnliuN - PruN)] 


/ 

m. 


m. 


27TmB 

\p< 


for Jb = \ 


— ) 3 ^^^ + P^n] for 


Q di 

jP -- 
'^^“2 


(26) 


where \pcm\ = (rnN+m.„'^y^g {rriN ujm = y^IPcmp + The isospin factors are 

In*^ttN = 3 and lA,A*^nN = 1 for the Ntt decay channels of N* and A* (or A) respectively. 
Putting the experimental values [2^ of Fb in Eq. (f2U|l , the values of coupling constants //have 
been fixed, which are shown in a Table (1). 


4 Results and discussion 

Now let us discuss the results of numerical evaluation. In Fig. ([2]), the imaginary (upper panel) and 
real (lower panel) part of pion self-energy for na (solid line) and vrp (dashed line) loops are plotted 

against the invariant mass M, where M = \J— k'^. The Landau and unitary regions for ttu and 
Trp loops are distinctly observed from the upper panel of the Fig. ([2]). The Landau regions of the 
TTfj and np loops are respectively (M = 0 to nia — m-,^ i.e. 0 to 0.25 GeV) and (M = 0 to nip — 
i.e. 0 to 0.63 GeV) whereas their unitary regions start respectively from (M = = 0.53 

GeV) and (M = nip + = 0.91 GeV) to oo. Since the unitary regions are far away from the 

pion pole (M = m.,^ = 0.14 GeV), the unitary cut contributions may not play any dominating 
role in the spectral profile of pion. Rather the Landau cut contributions seem to be responsible 
for determining the spectral profile of pion. For ttct loop, the broad vacuum spectral shape of a 
is taken into account by using Eq. (jl6p and the corresponding results are shown by dotted line in 
Eig. ([2]). Erom the upper panel of the Eig. ([2]), the dotted line indicates that due to the folding, 
the Landau and unitary cut regions are overlaped by each other and their contributions are also 
reduced. Eor the different set of temperature T and the three momentum k, the imaginary (upper 
panel) and real (lower panel) part of total self-energy after summing the contributions of ira and 
np loops have been shown in the Eig. ([3]). Here we see that the numerical strength of imaginary 
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Figure 4: Imaginary part of pion self-energy for different NB loops. B = 

A*(1620), iV*(1650), ^"*(1720) are displayed in lower panel while B = A(1232) is shown in upper 
panel, we have taken the medium parameters T = 0.15 GeV, = 0.453 GeV, for which the 
nucleon density become pN = 0.55/?o {po is saturation nuclear matter density). 



Figure 5: Same as Fig. (j3|) for rest of the baryons B = iV*(1440), iV*(1535), A*(1600) A^*(1700) 
(upper panel) and B = A*(1520), A*(1700) (lower panel). 
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Figure 6: The corresponding results of Fig. ([5]) for the real part of pion self-energy. Only NN loop 
contribution is additionally presented here which has vanishing Landau cut contribution in Fig. (11- 



Figure 7: The corresponding results of Fig. ([5]) for the real part of pion self-energy.. 
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Figure 8 : Imaginary (upper panel) and real (lower panel) part of total pion self-energy coming from 
the baryon loops are presented for different set of temperature T, baryon chemical potential /iw 
and three momentum k. 

part of pion self-energy is significantly enhanced in the low M regions (including the pion pole) for 
finite k. It also increases with the increasing of the k and T. The real part of the pion self-energy 
become very small around it’s pole position. 

Next we will discuss about the results of pion self-energy which are coming from the different 
NB loops. The Fig. (H)) demonstrates the imaginary part of pion self-energy for the baryons 
B = A(1232) (upper panel) and B = A*(1620), A*(1650), A*(1720) (lower panel). Whereas 
the results for baryons B = A*(1440), iV*(1535), A*(1600), iV*(1700) (upper panel) and B = 
A*(1520), A*(1700) (lower panel) are displayed in the Fig. ([5]). The result of the B = iV*(1710) 
is not revealed with the other baryons as its strength is very low. In those figures the Landau 
regions of different loops are clearly identified. As an example the Landau region of the NA loop 
is (M = 0 to m-A — mj\f i.e. 0 to 0.292 GeV). Fig. ([^ and ([7]) present the corresponding results of 
real part for different baryons. Summing all the baryon loops the total self-energy is represented 
in Fig. ([ 8 |) where upper and lower panel are exhibiting the imaginary and real part respectively for 
different set of {k, T and ^n)- Similar to mesonic loops, Imlls in low mass region is sensitively 
enhanced for hnite values of k and it also increases with the increasing of T and /tat. The upper 
panel of Fig. ([5]) clearly demonstrates this fact. The non-zero values of ImRsiM) for NB loops (or 
for 7 r 4 > loops) in the Landau cuts, 0 < M < ms — rriN (or 0 < M < m# — m-,^) basically 
interprets the statistical probability of pion due to scattering with medium constituents [27]. The 
possible scattering processes are as follows. During propagation in the medium, pion may disappear 
by absorbing a thermalized N (or tt) from the medium to create a thermalized B (or <I>). Again 
pion may appear by absorbing a thermalized B (or <h) from the medium as well as by emitting a 
thermalized N (or vr). Hence these kind of scattering processes provide the statistical probability 
to pion and this probability in the low M region is quite high for finite values of A:, T and //tv- 
From the upper panel of Fig. ([3|) and ([8|), one can notice that the low mass probability of pion for 
mesonic loops is dominating over that of the baryonic loops. 

To compare with some earlier results from Refs. miiaiie], we have generated the numerical 

values of pion optical potential , where ujk = \jk‘^ + rn‘^. Its imaginary (upper panel) 

and real (lower panel) part for mesonic (solid line) and baryonic (dotted line) matter are plotted 
against k in Fig. ([9|). For mesonic matter, our imaginary part of pion potential is qualitatively 
similar with the results, obtained by Shuryak [TT| (solid line with square) and Rapp et al. [15] 
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Figure 9: The imaginary (upper panel) and real part (lower panel) of pion potential for mesonic 
matter, obtained by Shuryak [T3] (solid line with square), Rapp et al. [15] (solid line with triangle) 
and for baryonic matter, obtained by Rapp et al. [16] (solid line with circle). Our corresponding 
results for meson loops, NA loop at A = 0.3 and 1.2 GeV are shown by solid, dashed and dash- 
dotted line respectively. 


(solid line with triangle) but quantitatively larger than them. One of the reason for this quantative 
difference may be that our results are originated from the effective nirp and tt-ttct interactions. In 
baryonic matter (considering N and A in the medium), our imaginary and real part both follow 
approximately similar k dependency as followed in Ref. [16| (solid line with circle). In Fig (|lUp . 
We have also compared our baryonic results with the results of Tolos et al. |18j . where the pion 
self-energy incorporates a momentum independent s-wave and a momentum dependent p-wave 
piece originated from the one and two nucleon-hole and delta-hole excitations, plus short-range 
correlations. Unlike to the results of Tolos et al. [18], our results and the results of Rapp et al. [T6] 
exhibit a non-monotonic momentum dependence of pion optical potential with a peak structure 
in its imaginary part. From the Eq. (|22p this nature can be grossly understood. The numerical 
integration of thermal distribution will roughly be close to a sine hyperbolic function of momentum, 
which will be a rapidly increasing function whereas the l/lGirk term before the integral will make 
it to be a decreasing function. Hence the net effect of these two opposite dependence of momentum 
create such kind of non-monotonic momentum dependent potential. In Ref. [18], incorporating 
the extra two loop kind of nucleon-hole and delta-hole excitations as well as an explicit function 
of A width with external variable may be one of the key reason for exhibiting the monotonically 
increasing behavior of pion optical potential with momentum. 

Using the imaginary and real part of pion self-energy in the Eq. ([6]), pion spectral function 
can be numerically estimated. The contributions of meson and baryon loops in the pion spectral 
function are cumulatively shown in Eig. m- After adding NA as well as the other NB loops 
with the meson loops, the peak of the pion spectral function is gradually reduced and shifted 
toward the higher M. This fact can be noticed very well in the lower panel of Eig. m for the 
parameters {T,fj,]\f,k) = (0.15,0.4,0) GeV. Erom the lower panel of Eig. ([8]), the positive values of 
ReH^ around the pion pole already indicates the fact. By changing the three momentum k from 0 
to 0.5 GeV, the peak strength of the pion spectral function has been reduced signihcantly, which 
can be observed in the upper panel of Eig. (ED- The reason for such kind of reduction is because 
of the enhancement of ImH at hnite k, which has been observed in the upper panel of Eig. ([3]) and 
([5]). This kind of softening of the pion pole for increasing momentum also observed in some earlier 
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Figure 10: The imaginary (upper panel) and real part (lower panel) of pion potential at temperature, 
T = 0.07 GeV and nucleon density, = Po- Solid (for A = 0.3 GeV) and dashed line (for A = 1.2 
GeV) show our results for A^A loop whereas the solid line with circle shows the corresponding 
results of Tolos et. al [18] . 


works mm- Prom Fig. CID, we also observe that the NA loop dominates over the rest of the 
NB loops to determine the pjy dependency of pion spectral profile. The Landau cut structures of 
the other NB loops are dominantly contributed to the higher M regions, which are slightly away 
from the pion pole. Hence, their contributions may be considered as a negligible correction part in 
the pion spectral function with respect to the contribution of iVA loop. Similar dominance of AA 
loop over the other baryon loops was also noticed by Post et al. |25j . However, by increasing the 
pion momentum k this correction part of ImH/m.^- can be dominant over the contribution of A A 
loop as shown in the Table (2). 


k 

ImH/mTr (AA) 

ImH/mTr {NB) 

ReH/m^ (NA) 

ReH/m^ (NB) 

0 

1.4 MeV 

1.4 MeV 

8.4 MeV 

11.4 MeV 

300 MeV 

38.6 MeV 

43.5 MeV 

15.1 MeV 

15.9 MeV 

500 MeV 

3.3 MeV 

10.1 MeV 

5.6 MeV 

3.4 MeV 


Table 2: Table shows imaginary and real part of pion self-energy coming from A A loop only and 
all NB loops (including AA) for different values of pion momentum k. 


Using the total pion self-energy after summing all the meson and baryon loops, we can get a 
full spectral function of vr meson at finite temperature and density. For different sets of (T, pjy) 
and k, this explicit structure of pion spectral function are respectively demonstrated in upper 
and lower panel of Fig. (|12|) . Due to noticeable Landau cut contribution of pion in the low mass 
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Figure 11: The contributions in pion spectral function from meson (dotted line), meson + A (dashed 
line) and meson + all baryon (solid line) loops are displayed for (T,fj,]\f,k) = (0.15,0.453,0) GeV 
(lower panel) and {T,fj,]\f,k) = (0.15,0.453,0.5) GeV (upper panel). 


region (below its pole position), there will be some Bose enhanced probability of p —)■ mr (along 
with TTTT —)■ p) below the vacuum threshold 2 m 7 r, which may have some contribution in the low 
mass dilepton enhancement. When we will apply this explicit in-medium pion spectral function 
to generate dilepton spectra via p meson modification, the probability of pion at high momentum 
will definitely be suppressed due to its statistical weight governed by the Bose-Einstein distribution 
function. Hence, a selective finite momentum range will be relevant in the dilepton spectra and the 
low mass enhanced probability of pion in that particular range may have some influence in the low 
mass dilepton spectra. This is our next interest of our future work. 

5 Summary and conclusion 

To summarize, the in-medium changes of pion propagation due to different mesonic and baryonic 
quantum fluctuation at finite temperature is investigated in the RTF. Since effective hadronic 
Lagrangian for ima and tttt/) interactions can alternatively describe the vacuum tttt scattering cross 
section, hence ttcj and irp loops have been considered as relevant mesonic quantum fluctuations of 
pion propagation in the effective theory. The imaginary and real part of pion self-energy for those 
mesonic loops have been evaluated in the RTF. To account the baryonic quantum fluctuation of 
pion propagation at finite temperature, an extensive set of nucleon-baryon loops are taken. After 
summing all the mesonic and baryonic loop contributions the total self-energy of pion in hot and 
dense nuclear matter has been numerically estimated. The imaginary and real part of this total 
self-energy determine the Breit-Wigner type structure of in-medium pion spectral function, which 
is appeared to be slightly complex in nature because of the non-trivial branch cut structures of self¬ 
energy function. At high temperature as well as density, the pion spectral profile are broadened with 
their attenuated peak structures. The main punch point of this article is the low mass probability 
of pion which is highly sensitive with pion’s momentum (k) and also with the medium parameters 
(r and pn)- This probability is completely appeared because of medium as Landau cuts of pion 
self-energy for different mesonic and baryonic loops are responsible for this contribution. Because 
of this low mass probability of pion, a non-zero Bose enhanced probability of p —)■ mr (along with 
Trvr —7- p) can be received below the vacuum threshold 2mT^, which promise to contribute in the low 
mass dilepton enhancement. However, a particular range of pion momentum will be responsible 
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Figure 12: Full pion spectral function for different sets of and k are shown in upper and 

lower panel respectively. 


for this low mass dilepton enhancement when we will apply this in-medium pion spectral function 
on the dilepton spectra via p meson modification. This efforts of calculating the low mass dilepton 
rate, which will be coming from the p meson self-energy for modified tttt loop, are in progress and 
will be addressed in due course. 

Acknowledgment : The work is financially supported by Fundacao de Amparo a Pesquisa do 
Estado de Sao Paulo, FAPESP (Brazilian agencies), under Contract No. 2012/16766-0. I am very 
grateful to Prof. Gastao Krein for his academic and non-academic support during my postdoctoral 
period in Brazil. 


References 

[1] R. Rapp, Adv. High Energy Phys. 2013, 148253 (2013); R. Rapp, J. Wambach, Adv. Nucl. 
Phys. 25, 1 (2000). 

[2] P. Mohanty, S. Ghosh, S. Mitra, Adv. High Energy Phys. 2013, 176578, (2013). 

[3] R. Arnaldi et al. for the NA60 collaboration, Phys. Rev. Lett. 100, 022302 (2008); Eur. Phys. 
J. C 61, 711 (2009); S. Damjanovic et al. for the NA60 Collaboration, J. Phys. G; Nucl. Part. 
Phys. 35, 104036 (2008). 

[4] R. Rapp, C. Gale, Phys. Rev. C 60, 024903 (1999) 

[5] S. Ghosh, S. Sarkar, S. Mallik, Eur. Phys. J. C 70 (2010) 251. 

[6] B. Friman and H.J. Pirner, Nucl.Phys. A 617 (1997) 496. 

[7] W. Peters, M. Post, H. Lenske, S. Leupold and U. Mosel, Nucl. Phys. A 632 (1998) 109. 

[8] S. Ghosh, S. Sarkar, Nucl. Phys. A 870 (2011) 94. 

[9] G. Chanfray and P. Schuck, Nucl. Phys. A 555 (1993) 329. 

[10] M. Herrmann, B. Friman and W. Norenberg, Nucl. Phys. A 560 (1993) 411 


15 






[11] M. Urban, M. Buballa, R. Rapp, J. Wambach, Nucl. Phys. A 673 (2000) 357. 

[12] A. Schenk, Nucl. Phys. B 363 (1991) 97. 

[13] G. Chanfray, D. Davesne, Z. Phys. A 349, 65 (1994). 

[14] E. V. Shuryak, Nucl. Phys. A 533 (1991) 761. 

[15] R. Rapp and J. Wambach, Phys. Lett. B 315 (1993) 220. 

[16] R. Rapp, J. Wambach, Nucl. Phys. A 573 (1994) 626. 

[17] R. Rapp, G. Ghanfray, J. Wambach, Nucl. Phys. A 617 (1997) 472. 

[18] L. Tolos, A. Ramos, A. Polls, Phys. Rev. C 65, 054907 (2002). 

[19] S. Mitra, S. Ghosh, S. Sarkar Phys. Rev. C 85 (2012) 064917. 

[20] A.J. Niemi, G.W. Semenoff, Ann. Phys. 152, 105 (1984) 

[21] R. L. Kobes and G. W. Semenoff, Nucl. Phys. 260, 714 (1985). 

[22] D. Fernandez-Fraile and A. Gomez Nicola, Fur. Phys. J. C 62, 37 (2009); Int. J. Mod. Phys. 
E 16, 3010 (2007); Eur. Phys. J. A 31, 848 (2007). 

[23] S. Ghosh, G. Krein, S. Sarkar, Phys. Rev. C 89, 045201 (2014). 

[24] S. Ghosh, Probing spectral properties of hadrons in hot and dense hadronic mat¬ 
ter (Ph.D. Thesis) HBNI (India), 2012, http://www.hbni.ac.in/phdthesis/thesis_june2013 
/PHYS04200704005_Sabyasachi_Ghosh.pdf 

[25] M. Post, S. Leupold, U. Mosel, Nucl. Phys. A 741, 81 (2004). 

[26] J. Beringer et al. {Particle Data Group) Phys. Rev. D 86, 010001 (2012). 

[27] H.A. Weldon, Phys. Rev. D 28, 2007 (1983). 

[28] P. Gonzalez, E. Oset, and J. Vijande , Phys. Rev. C 79, 025209 (2009). 

[29] S. Ghosh and S. Sarkar, Eur. Phys. J. A 49, 97 (2013). 


16 



